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ABSTRACT 

We present the evolutionary properties and luminosity functions of the radio sources belonging to the 
Chandra Deep Field South VLA survey, which reaches a flux density limit at 1.4 GHz of 43 /iJy at the 
field center and redshift ^ 5, and which includes the first radio-selected complete sample of radio-quiet 
active galactic nuclei (AGN). We use a new, comprehensive classification scheme based on radio, far- 
and near-IR, optical, and X-ray data to disentangle star-forming galaxies from AGN and radio-quiet 
from radio-loud AGN. We confirm our previous result that star-forming galaxies become dominant 
only below 0.1 mJy. The sub-millijansky radio sky turns out to be a complex mix of star-forming 
galaxies and radio-quiet AGN evolving at a similar, strong rate; non-evolving low-luminosity radio 
galaxies; and declining radio powerful (-P ^ 3 x 10^'* W Hz^^) AGN. Our results suggest that radio 
emission from radio-quiet AGN is closely related to star formation. The detection of compact, high 
brightness temperature cores in several nearby radio-quiet AGN can be explained by the co-existence 
of two components, one non-evolving and AGN-related and one evolving and star-formation-related. 
Radio-quiet AGN are an important class of sub-millijansky sources, accounting for ~ 30% of the 
sample and ~ 60% of all AGN, and outnumbering radio-loud AGN at < 0.1 mJy. This implies that 
future, large area sub-millijansky surveys, given the appropriate ancillary multi- wavelength data, have 
the potential of being able to assemble vast samples of radio-quiet AGN by-passing the problems of 
obscuration, which plague the optical and soft X-ray bands. 

Subject headings: galaxies: active — galaxies: starburst — radio continuum: galaxies — infrared 
radiation: galaxies — X-rays: galaxies 



1. INTRODUCTION 

The relationship between star formation and AGN in 
the Universe is one of the hottest topics of current extra- 
galactic research, at two different levels. On cosmological 
scales, the growth of supermassive black holes in AGN 
appears to be cor related with the growth of stellar mas s 
in galaxies fe.g.. IMerloni. Rudnick &: Di Matteoll2008[ ). 
On nuclear scales, the accreting gas feeding the black 
hole at the center of the AGN might trigger a starburst. 
The black hole, through winds and jets, can in turn feed 
energy back to its surroundings, which can compress the 
gas and therefore accelerate star formation but can also 
blow it away, thereby stopping accretion and star forma- 
tion altogether. The general consensus is that nuclear 
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activity plays a major role in the co-evolution of super- 
massive black holes and galaxies through the so-called 
"AGN Feedback" and indeed radio emission from AGN 
has been recently suggested to play an important role in 
galaxy evolution ^Croton et al. 2006). Moreover, radio 
observations afford a view of the Universe unaffected by 
the absorption, which plagues observations made at most 
other wavelengths, and therefore provide a vital contribu- 
tion to our understanding of this co-evolution. These two 
points imply that studies of the evolution of star- forming 
galaxies (SFG) and AGN in the radio band should pro- 
vide a better understanding of the link between the two 
phenomena. These are obviously done best by reaching 
relatively faint (< 1 mJy) flux densities, and hence the 
importance of characterizing the radio faint source pop- 
ulation. 
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After years of intense debate, the contribution to the 
sub-millijansky population from synchrotron emission re- 
suhing from relativistic plasma ejected from supernovae 
associated with massive star formation in galaxies ap- 
pears not to be overwhelming, at least down to ^ 50 /iJy, 
contrary to the (until recently) most accepted paradigm. 
Our deep (-S'i.4ghz > 43 /itJy) radio observations with the 
NRAO Very Large Array (VLA) of the Chandra Deep 
Field South (CDFS), complemented by a variety of data 
at other frequencies, imply in fact a roughly 50/50 split 
between SFG and AGN (Padovani et al. 2009), in broad 
agreement with other rec ent papers fe.g.. lSevmour et al.l 
[2008HSmorcic et al.l[2008l) . 

The purpose of this paper is to study the evolution 
and luminosity functions (LFs) of sub-millijansky radio 
sources through the VLA-CDFS sample. Apart from the 
topics mentioned above, this is important also for other 
issues, including: 

1. predictions for the source population at radio flux 
densities < 1 /xJy, which are relevant, for example, 
for the Square Kilometre Array (SKA). All existing 
estimates, in fact, had to rely, for obvious reasons, 
on extrapolations and are based on high flux den- 
sity samples. This affects particularly the highest 
redshifts, which can better be probed at fainter flux 
densities; 

2. the radio evolution of radio-quiet AGN. No radio- 
selected sample of radio-quiet AGN is currently 
available and this is badly needed to shed light on 
the mechanism behind their radio emission and al- 
low a proper comparison with radio-loud quasars; 

3. the fact that number counts by themselves do not 
necessarily reflect the relative intrinsic abundance 
of astrophysical sources, which req uires the deter- 
mina tion of the evolution and LF (jPadovani et al.l 
[2001 . 

We note that the evolution and LFs of sub-millijansky 
radio sources have been studied, so far, only in two fields; 
the COSMOS field (Smolcic et al. 2009a .bi and the Deep 
Spitzer Wide-area Infra Red Extragalactic (SWIRE) field 
(jStrazzullo et al.l I2010D . in both cases up to a maxi- 
mum redshift of 1.3 and without differentiating between 
radio-quiet and radio-loud AGN. Source classification 
was based on a rest-frame optical color scheme and on 
spectral energy distribution (SED) fitting to photometric 
data covering the UV to near-IR range respectively. 

We define as AGN sources in which most of the en- 
ergy is produced through physical processes other than 
the nuclear fusion that powers stars. In practice, this 
means that electromagnetic emission most likely related, 
directly or indirectly, to a supermassive black hole is pre- 
dominant in at least one band. A small fraction of AGN 
have, for the same optical power, radio powers three to 
four orders of magnitude higher than the rest. These 
are called "radio-loud" quasars and most of the energy 
they emit is non-thermal and is associated with powerful 
relativistic jets, although thermal components associated 
with an accretion disk may also be observed, especially 
in the optical/UV band. Radio galaxies are also char- 
acterized by strong radio jets (manifested also through 
radio lobes), typically laying in or near the plane of the 



sky, and a fraction of them (the most powerful ones) 
are thought to be radio-loud quasars, which have instead 
their jets oriented with a small angle to the line of sight 
(e.g.. lUrrv fc Padovani|[T9 95). We define as "radio-quiet" 
AGN in which jets are either not present or make a tiny 
contribution to the total energy budget over the whole 
electromagnetic spectrum, which is dominated by ther- 
mal emission. All other AGN we call radio-loud. Note 
that radio-quiet AGN are not radio-silent. Indeed, the 
radio power of many low-luminosity radio galaxies, the 
so-called Fanaroff- Riley (FR) Is ("low-power radio- loud 
AGN" according to our nomenclature) overlap with that 
of radio-quiet AGN, which can generate some confusion 
and requires great care during the classification process, 
which needs to involve also the X-ray and far-IR bands 
(Section 12. 4p . However, the two classes are physically 
distinct (see Sections 15. 31 and [5. 8p . although the origin of 
radio emission in radio-quiet AGN is still not clear (but 
see Section [STT]) . 

Translating these high-level definitions into a classi- 
fication scheme requires a wealth of multi-wavelength 
data, which wer e described in our previous papers. 
IKellermann et al"] (|2008) (Paper I) presented the radio 
data of the VLA-CDFS sample, to gether with optical im - 
ages and X-ray counterparts, while lMainieri et all (|2008[ ) 
(Paper II), discussed the optical and near IR counter- 
parts to the observed radio sources and, based on rest- 
frame colors and the morphology of the host galaxies, 
found evidences for a change in the sub-millijans ky ra- 
dio source population below « 80 microJy. Tozzi et al.l 
(200 3) (Paper I I I), dea lt with the X-ray properties, while 
I Padovani et al.l (|2009[ ) (Paper IV), discussed the source 
population. This turned out to be made up of SFG and 
AGN at roughly equal levels, with the AGN including ra- 
dio galaxies, mostly low-power (FR Is), and a significant 
{^ 50%) radio-quiet component. Paper IV made also 
clear that the "standard" definitions of radio-loudness, 
based on radio-to-optical flux density ratios, R, and ra- 
dio powers, were insufficient to identify radio-quiet AGN 
when dealing with a sample, which included also star- 
forming and radio-galaxies, as both classes are or can be 
(respectively) characterized by low R and radio powers as 
well. R, for example, is useful for quasar samples, where 
it can be assumed that the optical flux is related to the 
accretion disk, but obviously loses its meaning as an in- 
dicator of jet strength if both the radio and the optical 
band are dominated by jet emiss ion, as might be the case 
in FR Is (jChiaberge et al.lll999f) . Source classification in 
Paper IV was then based on radio, optical, and X-ray 
data, and was meant to provide a robust upper limit to 
the fraction of SFG at sub-millijansky levels. SFG can- 
didates were selected based on their (low values of) R, 
(low) radio power, (non-elliptical or SO) optical morphol- 
ogy, and (low) X-ray power (ix)- The fact that X-ray 
upper limits above the AGN threshold (10''^ ergs s~^) 
were also included was conservative in the sense that it 
maximized the number of SFG, as some of these sources 
could still be AGN. Furthermore, the selection of radio- 
quiet AGN candidates was only approximate, as it was 
based solely on R and Lx and suffered from uncertain- 
ties in the optical K-correction, possible contamination 
by radio galaxies, and the exclusion of X-ray upper lim- 
its (see Paper IV for details). In order to deal with the 
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Fig. 1. 

formation. The solid line is the best fit, while the dashed lines 
represent the scatter. See text for details. 

evolution and luminosity functions of the various classes 
of sources, we need to refine our classification. In partic- 
ular, the CDFS field has been observed by Spitzer and 
therefore near- fSection l2.2p and mid/far-IR (Section l2.3l) 
data are available for our sample. The source classifica- 
tion used in this paper relies then on a combination of 
radio, IR, optical, and X-ray data (Section [Z4|) . 

Section 2 describes the updated classification of the 
VLA-CDFS sample, while Section 3 studies its evolu- 
tion. Section 4 derives the LFs for various classes, while 
Section 5 discusses our results. Finally, Section 6 sum- 
marizes our conclusions. Throughout this paper spectral 
indices are defined by 5*^ ex i'^°' and the values Hq = 70 
km s^^ Mpc^^, 51m — 0.3, and JIa = 0.7 have been used. 

2. THE UPDATED SAMPLE CLASSIFICATION 

2.1. Redshifts 

Our sample includes all VLA-CDFS sources with reli- 
able optical counterparts and eight empty fields, for to- 
tal of 256 objects, 193 of which belong to a complete 
sampleQ (see Paper IV for details). 92% of the sources 
in the complete sample now have redshift information 
(74% spectroscopic) as compared to only 77% in Pa- 
per IV. We have included in this paper spectroscopic 
redshifts from a dedicated follow-up program performed 
with the VIMOS spectrograph at the VLT (Bonzini et 
al., in preparation). We also used recently published red- 
shifts for the co unterpar ts of Chandra sourc es in this field 
(jTreister et al.r 2009: Sil verman et al.ll2010D and the pho- 

^ Four more sources belonging to the complete sample have very 
uncertain counterparts (see Paper II) and for one other source, 
very close to a bright star, we could not get reliable photometry. 
The inclusion of these sources in any of the classes described below 
would change our results by much less than la. 



tometric redshifts published by the Multiwavelength Sur - 
vey by Yale-Chile (MUSYC) (ICardamone et al.llMol) . 
which are based on 32 photometric bands. 

As shown in Fig. [1] redshift is strongly correlated 
with magnitude, albeit with some scatter. The best 
and simplest approach to estimate the redshift for the 
16 objects in the complete sample without observed red- 
shifts is then to derive it from their magnitude by using 
the relationship shown in the figure (solid line), that is 
logz = 0.166Vtnag — 3.85. This was derived applying 
to the w hole sample the o rdinary least-square bisector 
method (jlsobe et al.l 119901 ) , which treats the variables 
symmetrically. Including only spectroscopic redshifts, 
or only the complete sample, or excluding sources with 
likely AGN contamination in the optical band (based on 
iSzokolvet al. 2004), all give relations within la from the 
adopted one. The effect of this assumption on our results 
is discussed in Section 15.11 

Note that while objects well to the left of the corre- 
lation can be explained as having an AGN component 
in the optical band, the single source in the lower right 
part of the diagram is ~ 7 magnitudes fainter than the 
average and therefore well into the dwarf galaxy regime. 
However, its photometry is affected by its closeness to 
a bright star, which might explain at least in part its 
faintness. 

2.2. Near-IR data 

The usage of S yitzer Infrared Array Camera (IRAC; 
iFazio et all |2004| ) colors to identify AGN has been 
discu ss ed at length in the l i teratu r e fe.g.. iLacv et al.l 
20041: iHatz iminaoglou et all [20051: iStern et al.l 120051: 
Saiina et a l. 2005; Card amone et al.|[2008D . Although it 



is by now evident that only some classes of extragalac- 
tic sources occupy restricted regions of parameter space 
in such plots, it is nevertheless also clear that there are 
broad trends which can be used to, for example, identify 
possible misclassifications. 

Fig. [2] plots the IRAC flux density ratios S'8.o/'S'4.5 
vs. 55.8/5*3.6 for our sources classifled as in Paper IV, 
where the fiux densities refer to all four IRAC chan- 
nels at 3.6, 4.5, 5.8, and 8.0 /im. The IRAC data come 
from the Spitzer IRAC /MUSY C Public Legacy surv ey in 
the Extended-CDFS (SIMPLE: iDamen et al.ll20lll) . We 
cross-correlated the SIMPLE catalogue with the VLA- 
CDFS sources, accepting matches with separations less 
than 2". The SIMPLE catalogue has convolved the im- 
ages associated with each IRAC channel to match that 
of channel 4 (8.0 /^m), the one with the lowest resolu- 
tion, so that reasonably accurate colors may be obtained 
from the four IRAC bands. We have used total fluxes 
and applied the prescribed normalization to produce flux 
densities in /iJy. Fig. [2] shows the following: (a) most 
AGN candidates fall around the locus of sources whose 
mid-IR spectrum can be characterized by a single power 
law (dotted line); a significant number of AGN is also 
within the so-called "Lacy 's wedge" (dashed lines), which 
is where most unobscured, broad-lined (type 1) ACN are 
thought to lie (Lacy et al. 2004) . Note that highly ob- 
scured sources r night also occupy that region of pa ram- 
eter space (e.g.. lDasvra et al.ll2009HPrandonill2010D : (b) 
most SFG candidates are distributed in a vertical band 
centered around S^.s/S^s '-^ 0.6 — 0.8, which is where 
polycyclic aromatic hydrocarbon (PAH)- and starlight- 
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Fig. 2. — IRAC color-color plot for the AGN and SFG candidates 
selected in Paper IV and four FR Is from the SWIRE field. The 
dotted line indicates the locus of sources whose spectrum can be 
described as a power-law over the four IRAC bands. The dashed 
lines indicate the so-called "Lacy's wedge" , which is where most 
AGN are thought to lie. The solid lines denote a more restrictive 
region, which takes into account the fact that for z > 0.5 PAH- 
and star-ligh t-dominated sources can be inside "Lacy's wedge" 
IIDasyra et a l. 2009). See text for details. 



domi nated sources are expected to lie fe.g.. lSaiina et alJ 
[2005h : (c) we also plot four "bona fid e" FR I from the 
SWIRE field (|Vardoulaki et al. 1 120081) . which fall in the 
region where galaxies with an ol d stellar population are 
located (e.g., iSaiina et al.l 120051 ) . We have quite a few 
sources in the same area, which is consistent with one of 
the main results of Paper IV, that is the dominance of 
low-luminosity radio galaxies amongst radio-loud AGN. 
This is reassuring and shows that the SFG /AGN di- 
vision derived in Paper IV is overall correct. The most 
interesting features in Fig. [2l however, are the excep- 
tions to the above, namely: (a) the eight AGN candi- 
dates in the top left part of the diagram; these are all 
but one at low redshift (z < 0.25), low radio power 
(logPi.4GH^ < 22.6), low X-ray (2 - 10 keV) power 
(logLx < 41.6) sources, which had been classified as 
AGN solely because their optical morphology was SO (5) 
or elliptical (2). A closer look at their images shows 
that two of them show (weak) signs of spiral arms and 
four more (all SO) have only low-resolution Wide Field 
Imager (WFI) data, which means that the presence of 
spiral arms cannot be e xcluded. Their loc ation in the 
PAH-dominated region (jSaiina et al.l I2005D suggests a 
re-classification as SFG for all of them apart from one 
AGN with two spiral galaxies at a distance of ~ 3", 
which means its IRAC fiux is most likely contaminated 
(its rest-frame radio-to-optical fiux density ratio is also 
~ 2, which is typical of radio-loud AGN: see below); 
(b) the ten SFG candidates with S^.s/S^s > 1 and 
1 < 5'8.o/'S'4.5 < 3; this is more restrictive than the 



Lacy's wedge as it takes into account the fact that for 
z > 0.5 PAH- and star-light-dominate d sources can be 
inside that wedge ()Dasvra et al.l I2OO90 . Most of these 
sources have X-ray upper limits larger than 10''^ erg/s, 
which makes sense since this was one of the reasons they 
were classified as SFG in the first place. The location 
of these sources suggests a re-classification as AGN. In 
summary, seven sources were re-classified from AGN to 
SFG and ten sources previously classified SFG are now 
classified as AGN. 

2.3. Far-IR data 

It is well know that the global far-IR and radio 
emission are tightly a nd linearly correlat ed in star- 
forming systems (e.g., I Sargent et al.l 120101 and refer- 
ences therein). This is usually expressed through the 
so-called q parameter, that is the logarithm of the 
ratio of far-IR to radio flux density, as defined by 
iHelou. Soifer. fc Rowan-RobinsonI ()1985D . We take ad- 
vantage of the relatively narrow dispersion of q for star- 
forming systems to further refine our SFG/AGN separa- 
tion and also to improve on our radio-quiet - radio-loud 
AGN division, as the latter do not follow the IR - radio 
correlation typical o f SFG (e.g-. lSopp fc AlexandedlTQQTI : 
iSargent et al.ll20Tol) . This is vital to separate radio-quiet 
AGN from radio galaxies, as R is not very useful in this 
case (Section[T]) and, like radio-quiet AGN, radio galaxies 
can also have relatively large X-ray powers. 

We have used a catalog of 70 /xm Multiband Imaging 
Photometer for Spitzer (MIPS) flux densities from the 
Far- Infrared Deep Extragalactic Survey (FIDEL; Dick- 
inson et al. in preparation) for our evaluation of q. We 
cross correlated the VLA-CDFS radio sources with the 
FIDEL catalogue using a radius of 8" (about half the 
Spitzer 70 fim point-spread-function). For those sources 
undetected by the FIDEL survey (but still within the FI- 
DEL coverage), we assume an upper limit of 2.5 mJy as 
this is approximately the 5cr survey limit. To these data 
we add 24 /im flux densities from the Great Observato- 
ries Origins Deep Survey (GOODS) whenever available, 
and thus we obtain SEDs sampled at up to eight wave- 
lengths: 20 and 6 cm in the radio from our VLA surveys; 
70 A*m and 24 ^m in the IR from FIDEL and GOODS; 
and 8.0 /im, 5.6 fim, 4.5 /im, and 3.6 /im in the near-IR 
from SIMPLE. 

We then pro c eeded to find the template SED from 
the iDale et al.l (|2001l ) SFG models that best matches 
the Spitzer data. We use the source redshifts to place 
each of the 64 models into the observed frame for that 
source, and set the normalization by requiring that each 
model SED pass through the measured 70 /im flux den- 
sity for that galaxy. This, in effect, places an extra 
weight on 70 fim data since it is our only measure- 
ment of the smooth modified blackbody portion of the 
SED. We then select the model which minimizes the least 
squares fit to the photometry of the four IRAC chan- 
nels and the MIPS 24 /im data (when available). Once 
the best-fitting model has been selected, we derive the 
rest-frame 60 /im and 100 /xm flux densities to determine 
FIR, the total far-IR flux between 42.5 /i m and 122.5 /im 
(|Helou. Soifer. fc Rowan-Robinsonlll985| ) 



FIR = 1.26 X 10-^* [2.58/6o^„ -1- f, 



lOO^im. 



W m 



(1) 
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where the flux densities, /, are in Jy. Similarly, we con- 
vert the observed 1.4 GHz radio emission to the rest 
frame using its measured spectral index between 1.4 GHz 
and 4.86 GHz, where available (^ 80% of the sample: see 
Paper I), or assuming a spectral index a^ = 0.7 (the mean 
of the sample) otherwise. The value of q is then calcu- 
lated as the logarithm of the ratio of far-IR to 1.4 GHz 
flux density: 



q = log[{FIR/3.75 x 10^^)/Siaghz] 



(2) 



where the numeric factor is the frequency in Hz corre- 
sponding to a wavelength of 80 /xm. 

Given the large fraction (^ 50%) of upper limits on 
q, one cannot readily look for a bimodality in its dis- 
tribution to define a dividing line between star-forming 
and non star-forming sources. The median of the detec- 
tions should be however quite well defined, as its value 
is 2.16 and most upper limits are below 2.2. Since 96% 
of the detections above the median are below 2.64 and 
assuming a symmetric distribution, one finds a lower 
end at around 2.16 - (2.64 - 2.16) - 1.7. We then 
assume in the following that sources characterized by 
g > 1.7 are star-formerqj (upper limits above this value 
excluded). This is the same dividing value assumed by 
iMachalski fc CondonI ()1999() . Twenty- two of our candi- 
date SFG have q < 1.7 and therefore cannot be star- 
forming systems. These were then re-classified as radio- 
loud AGN. These sources fall in the region where passive 
galaxies are found in the IRAC color-color plot, which is 
consistent with this re-classification, given that most of 
our radio-loud AGN should be radio galaxies. 

Finally, eight radio-quiet AGN candidates were found 
to have q < 1.7, while nineteen radio-loud ones had 
q > 1-7, which reflects the approximation of our pre- 
vious classification. These objects were re-classified as 
radio-loud and radio-quiet respectively. 

2.4. Revised classification 

To summarize, based on the results presented in Pa- 
per IV and in the previous sub-sections, our candidate 
star-forming galaxies are defined as fulfilling the follow- 
ing initial requirements: 

1. R ^ log(5i.4GHz/S'v) < 1.7 (where Sy is the V- 
band fiux density) 

2. P, < 1024-5 W Hz-i 

3. optical morphology different from elliptical or 
lenticular 

4. Lx(2-10 keV) < 10"^ ergs s^^ for X-ray detections, 
no limit otherwise. 

As discussed in Paper IV, the first two criteria include 
^ 90% of spirals and irregulars, the third one excludes 
sources not associated with star formation at our red- 
shifts ((z) ^ 1.1), while the fourth one excludes AGN. 
These are then supplemented by the following additional 
requirements, which can overrule the previous ones if nec- 
essary: 

^ Our results are only weakly dependent on this choice. For 
example, if we defined as star-formers sources with q > 1.8 our 
SFG complete sample would only lose three objects, a 4% effect 
(see Tab. tg. 
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A flow chart of our classification scheme. See text for 



5. IRAC constraints: the region of parameter space 
defined by S's.s/^'s.e > 1 and 1 < S'8.o/S'4.5 ^ 
3 {S5.s/S3.6)°'^^ (AGN region) is excluded; sources 
not classified as SFG by the previous criteria but 
with 0.45 < S^.a/Sa.e < 1-0 and S'8.0/5'4.5 > 2.5 
(PAH-dominated region) are also included (Sec- 
tion [2^ 



6. MIPS constraints: q > 1.7, upper limits above this 
value excluded f Section [23t 



Note that constraints number 2 and 3 have become 
almost irrelevant for our classification given these two 
new requirements. Objects not fulfilling this sequence of 
criteria are considered to be AGN. Radio-quiet AGN are 
defined initially as follows: 

1. R< 1.4 

2. Lx(2 - 10 keV) > lO'*^ ergs s^^ (detections only) 

As discussed in Paper IV, the first criterion is the "clas- 
sical" definition of radio-quiet AGN converted to the 1.4 
GHz and V bands. These are then supplemented by the 
following additional requirements, which can overrule the 
previous ones if necessary: 

3. IRAC constraints: the region of parameter space 
defined by S'8.o/'S'4.5 > 2.5 and 0.45 < S'5. 8/5*3.6 < 
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1.0 (PAH-dominated region) is excluded (Sec- 
tion [2^ 



4. MIPS constraints: q > 1.7, upper limits above this 
value excluded (Section 12.31) 

AGN not fulfilling this sequence of criteria are consid- 
ered to be radio-loud. Fig. [3] summarizes our classifica- 
tion scheme. We note that, while it is relatively simple to 
distinguish SFG from radio-loud AGN due to their differ- 
ent q values, the situation is more complicated when one 
has to differentiate SFG from radio-quiet AGN. This is 
done also based on the location on the IRAC color-color 
plot, which highlights the obvious outliers, but mostly 
on the basis of X-ray power. Since many SFG still have 
upper limits on Lx(2 — 10 keV) > 10^^ ergs s~^, we can- 
not exclude that more radio-quiet AGN are present in 
our sample, especially amongst the SFG with the high- 
est limits on Lx(2 — 10 keV), which tend to be at higher 
redshifts. The inclusion of deeper X-ray data in our anal- 
ysis will help us sort out this issue. 

It is instructive to see how representative local sources 
get classified by our scheme. To this aim, we have 
used the NASA/IPAC Extragalactic Database (NED) 
and NASA's Astrophysics Data System (ADS) to get 
multi-wavelength data for a few objects. For example, 
NGC 1068, the prototype Seyfert 2 galaxy, residing in a 
spiral host, has low R, P^, and Lx(2 — 10 keV) values. 
Coupled with a g > 1.7 it would then be classified as an 
SFG but its location on the IRAC color-color plot puts 
it firmly with the radio-quiet AGN. NGC 1052, an el- 
liptical galaxy often classified as a low-ionization nuclear 
emission- line region (LINER), has also low R, P^-, and 
Lx(2 — 10 keV) values but its low q makes it a radio-loud 
AGN. M 82, the prototype starburst galaxy, again has 
low i?, Pr, and Lx(2 - 10 keV) values but its q > 1.7 
and location on the IRAC diagram classify it as a SFG. 
And NGC 1275, a cD (radio) galaxy at the centre of the 
Perseus cluster, which in the literature has been classi- 
fied, amongst other things, as a Seyfert 1.5 and a blazar, 
with its high R, Pj., and very low q is undoubtedly a 
radio-loud AGN. 

Table [1] and Fig. |4] present the Euclidean normalized 
number counts for the revised sample, compared to those 
from Paper IV shown in the figure. As expected from the 
revised classification, SFG show a small decrease, while 
AGN and radio-quiet AGN increase slightly in number. 
However, the revised number counts are still within la 
from the old ones, and most of the largest changes hap- 
pen at higher flux densities. Note that SFG are still pre- 
dominant below ~ 0.1 mjy, which is also the flux density 
at which radio-quiet AGN start to outnumber radio-loud 
ones. 

AGN make up 49^1%(where the errors are based 
on binomial statistics: iGehrelsl 119861 ) of sub-millijansky 
sources and their counts are seen to drop at lower flux 
densities, going from 100% of the total at ^ 10 mJy 
down to 41% at the survey limit. SFG, on the other 
hand, which represent 50_y% of the sample, are missing 
at high flux densities but become the dominant popula- 
tion below w 0.1 mJy, reaching 59% at the survey limit. 
Radio-quiet AGN represent 28^5% (or 57% of all AGN) 
of sub-millijansky sources but their fraction appears to 
increase at lower flux densities, where they make up 84% 
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Fig. 4. — The Euclidean normalized 1.4 GHz CDFS source 
counts: total counts (black triangles), SFG (filled green circles), 
all AGN (red squares), radio-quiet AGN (open blue circles), and 
radio-loud AGN (open magenta triangles). Error bars correspond 
to la errors (Gehrols 1986). The dotted symbols are the source 
counts from Paper IV (shifted by 0.035 dex for clarity). The 
long-dashed and dot-dash ed lines are t he ra dio-q uiet AGN num- 
ber co unts predicted by IWilman et al.l I I200SI) and I Padovani et al.l 
II2009I ) respectively, based on X-ray data. See text for details. 

of all AGN and ~ 34% of all sources at the survey limit, 
u p from fti 5% at » 1 mJy . 

iMiddelberg et aH (j2011t ) have recently detected with 
the Very Long Baseline Array (VLBA) 20 VLA-CDFS 
sources using a resolution of ~ 0.025". With a limit 
of ^ 0.5 mJy, very long baseline interferometry (VLBI) 
detections above z > 0.1 are most likely to be due to 
AGN. Reassuringly, aU of the 20 detected VLBA objects 
(which have z > 0.15) were classified as radio- loud AGN 
by our method. 

Fig. [5] shows radio power versus redshift for our 
sources, with the dotted lines indicating 43 /LtJy, the 
faintest radio flux density of our sample (lower line) , and 
100 fiJy (upper line: see Section [3]) for a^ = 0.7. 

3. EVOLUTION 
3.1. Ve/Vg, analysis 

We first study the evolutionary properties of the VLA- 
CDFS sa mple through a variati on of the WKnax test 
(ISchmidtJ[l968l). td ie K/K test (lAvni fc Bahcal]l " fT980l : 
iMorris et al.lll99ll) , that is the ratio between enclosed 
and available volume. This is because we do not have 
a single flux limit but the sensitivity of our sample is a 
function of the position in the held of view (see Paper 
I) and, consequently, the area of the sky covered at any 
given flux density (usually known as the sky coverage) 
is flux dependent. This ranges from a maximum of 0.2 
deg^ for radio flux densities > 295 mJy, to 0.14 deg^ at 
100 mJy, down to 0.01 deg^ at the flux density limit. 

Values of (14/14) significantly different from 0.5 and a 
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Fig. 5. — Radio power versus redshift for our sample. Filled 
circles indicate SFG, filled squares radio-loud AGN, and open cir- 
cles radio-quiet AGN. Crosses denote redshifts estimated from the 
optical magnitude. The power - redshift relationships for a fiux 
density of 43 fijj (lower dotted line) and 100 /ijy (upper dotted 
line) assuming a radio spectral index of 0.7 are also shown. 

distribution significantly different from uniform indicate 
evolution, which is positive (sources were more luminous 
and/or more numerous in the past) when (T4/V^) > 0.5, 
and negative (sources were less luminous and/or less nu- 
merous in the past) when < 0.5. Moreover, one can 
also fit an evolutionary model to the sample by finding 
the evolutionary parameter which makes (Vc/Vg) = 0.5. 
Note that the T4/T4 test is independent of the shape of 
the luminosity function, unlike the maximum likelihood 
method used below. 

We have computed Ve/Va, values for our sources tak- 
ing into acc ount the appropriate s ky coverage (see eqs. 
42 and 43 of lAvni fc Bahcal l' 1980), k-correcting the ra- 
dio powers as described in Section 12.31 Statistical er- 
rors are given by cr == 1/V12 N (jAvni k BahcalllllQSOO . 
We estimate the significance of the deviation from the 
non-evolutionary case by assessing the probability Pev 



that the T4/T4 distribution is different from uniform ac- 
cording to a Kolmogorov-Smirnov (KS) test. Similar 
results are obtained by using the deviation from 0.5 of 
(Ve/Va). To have a simple estimate of the sample evo- 
lution we have also derived the best fit parameter k^ 
assuming a pure luminosity evolution (PLE) of the type 
P{z) = P(0)(1 -I- z)''^ or a pure density evolution (PDE) 
of the type $(z) = $(0)(1 -|- z)''", where <i>(z) is the 
luminosity function. 

We assume that some luminosity evolution takes place, 
based on previous studies in the radio and other bands. 
When the best fit indicates negative luminosity evolu- 
tion (i.e., fci < 0), however, we fit a pure density evolu- 
tion model as well, which we feel is more physical in this 
case. Note that for a single power law LF $(P) oc P~^ 
the evolutionary parameters in the two cases are related 
through the simple relationship kjj = ^^(7 — 1) (e.g., 
IMarshanetHirTosl . 

Our results are shown in Table [21 which gives the 
sample in column (1), the number of sources in column 
(2), the mean redshift in column (3), the percentage of 
sources with redshift estimated from the magnitude in 
column (4), (14/14) in column (5), the probability Pev 
that the 14/T4 distribution is different from uniform in 
column (6), and the best fit parameters k^ and kjj (when 
applicable) in column (7) and (8) respectively (only when 
Pev > 95%). The mean redshift is calculated taking into 
account the effect of the sky coverage, that is each ob- 
ject is weighted by the inverse of the area accessi ble at 
the fl ux density of the source (see, e.g., iPadovani et alJ 

The fractional redshift distributions for the different 
classes are shown in Fig. [HI As also shown in Table 
[2I radio-quiet AGN have the highest (z) and reach z ~ 
3.9. Radio-loud AGN have the broadest distribution, 
extending up to 2; ~ 6 but with a lower (z) than radio- 
quiet ones. SFG, despite the relatively narrow redshift 
range (z < 2.3) have a (z) not too different from that 
of radio-loud AGN and much smaller than radio-quiet 
ones. Most classes show strong peaks at z ^ 0.5 — 0.75. 
These are related for all classes apart from radio-quiet 
AGN to t he two larg e-scale structures detected in the 
CDFS by IGilh et al.i (|200l in the 0.664 < z < 0.685 
and 0.725 < z < 0.742 ranges, which, as shown in the 
figure, contribute substantially to the observed peaks. 
The effect of these structures on our results is discussed 
in Section [Ol 

The main results on the sample evolution are the fol- 
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TABLE 2 

Sample Evolutionary Properties: Vo/Vk analysis 



o 
o 

CO 



Sample 




N 


(z) 


estimated z % 


(VJV,) 


^GV 


fcL" 


fe^b 


All sources 




193 


1.18 ±0.07 


8.3% 


0.544 ±0.021 


98.8% 


n q+O-3 

u.y_o.4 




Star-forming galaxies 




71 


0.90 ±0.07 


1.4% 


0.655 ±0.034 


> 99.9% 


9 c+0.2 
^•■^ 0.3 




All AGN 




122 


1.44 ±0.09 


12.3% 


0.479 ± 0.026 


65.6% 


c 


c 


Radio-quiet AGN 




36 


1.73 ±0.16 


2.8% 


0.727 ± 0.048 


> 99.9% 


o c:+0.2 




Radio-loud AGN 




86 


1.26 ±0.11 


16.3% 


0.375 ±0.031 


99.5% 


-s.otli 


-1.6 ±0.4 


Radio-loud AGN, P < lO^'*'^ 


WHz-i 


53 


0.84 ±0.06 


9.4% 


0.432 ± 0.040 


76.2% 


c 


c 


Radio-loud AGN, P > lO^^-S 


WHz-i 


33 


2.01 ±0.26 


27.3% 


0.285 ±0.050 


99.9% 




-^■'>-0.9 



Pure luminosity evolution P(z) = P(0)(1 -I- z)'''^ 
' Pure density evolution $(z) = $(0)(1 + z)''o 
Pcv < 95%: no evolution required 
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redshift 

Fig. 6. — Fractional redshift distributions for the different classes 
of sources, deconvolved with the appropriate sky coverage. The 
dashed areas denote redshifts estimated from the optical magni- 
tude. Error bars represent the Itr range based on Poisson statistics. 
The dotted lines exclude all sources in the two large-sc ale concen- 
tratio ns at 0.664 < z < 0.685 and 0.725 < z < 0.742 IIGilli et al.l 
[2003) . 



lowing: 



1. The whole sample has (K/Ki) > 0.5 and shows 
a significant departure from the non-evolutionary 
case {Pcv ~ 99%) with an evolution characterized 
by kL = 0.9t"ol 

2. SFG evolve at a very high significance level (Pcv > 
99.9%); their evolutionary parameter for the case 



of pure luminosity evolution is k^ =2.5 



+0.2 
-0.3- 



AGN as a whole do not appear to evolve, as their 
(Ve/Vs.) is only slightly below 0.5 (by ~ 0.8a) and 
the Vc/Vg, distribution is not significantly different 
from uniform (Pev ~ 66%). 



4. Radio-quiet AGN, however, evolve very signifi- 
cantly (Pcv > 99.9%) with kL = 2.5loJ, the same 



value as that of SFG. 



1.3' 



5. Radio- loud AGN also evolve significantly (Pcv = 
99.5%) but in the negative sense, with fc^ — 
-3.0tJ;? or ko = -1-6 ± 0.4. However, this is 
largely due to the high power sources. Fig. [7] 
shows that T4/T4 values are strongly dependent 
on radio power, with (T4/T4) becoming signifi- 
cantly (- 2.8a) smaller than 0.5 for P > 2 x lO^" 
W Hz^^. We then split the radio-loud AGN sam- 
ple at P = 1024-5 W Hz-i. The low-power sub- 
sample does not appear to evolve, as its (K/Va) is 
not significantly (^ 1.7cr) < 0.5 and its Vc/Va dis- 
tribution is not significantly different from uniform 
(Pev ^ 76%). On the other hand, the high-power 
sub-sample anti-evolves at a very high significance 
level (Pev ~ 99.9%), with ko = -l-Sto.g- Because 
of the luminosity difference, the two sub-samples 
have also very different redshift distributions, with 
(z) - 0.8 (range: - 0.1 - 2.3) and 2.0 (range: 
^ 0.5 — 5.8) respectively. This also means that the 
low-redshift evolution of the high-power sub-class 
is not well determined as, for example, only eight 
such sources (24%) have z < 1 (see Section |4]) . 



3.2. Maximum Likelihood Analysis 

A more general approach to estimate the evolution, 
and at the same time the LF, is to perform a max- 
imum likelihood fit of an evolving luminosity function 
to the observed distribution in luminosity and redshift. 
This approach makes maximal use of the data and is 
free from arbitrary binning; however, unlike the 14/14, 
test, it is mode l depen dent. We follow the prescription of 
[Marshall et al.l ()1983[ ) and minimize the following quan- 
tity 



N 



^ = -2^1n[$(P„z,)]-h2 



<i>{p,z)n{f{p,z))^dzdP 

az 

(3) 



where $(P, 2;) is the luminosity function, fi(/) is the 
sky coverage, and dV is the differential comoving vol- 
ume. The sum is extended over the whole sample. 
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TABLE 3 

Sample Luminosity Functions and Evolution: 



MAXIMUM LIKELIHOOD ANALYSIS 
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Model 


71 
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logP, 
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Star-forming galaxies 
All AGN 




PLE 
PLE 
PLE 
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Fig. 7. — Vc/Vk values and means versus radio power for radio- 
loud AGN. The dashed line indicates the mean value for the whole 
sample, while the dotted line is the non-evolutionary value (0.5). 
Vertical error bars represent the la range based on Poisson statis- 
tics while the horizontal ones give the covered range of powers. 

while the double integral is computed over the luminos- 
ity range appropriate for the ad opted evolution and o ver 
the observed redshift range (see lMarshall et al1ll983l for 
more details). The best fit parameters are determined 
by minimizing S and their associated errors are com- 
puted by varying the parameter of interest until an in- 
crement AS* over the minimum value is obtained, la 
errors for one parameter correspond to AS* =1.0 "while 
confidence contours for 1, 2, and 3a levels for two in- 
teresting parameters are derived for AS" = 2.3, 6.17, 
and 11.8 respectively (jPress et al.| [l986V We consider 
one and two power-law LFs, that is $(P) oc P"7-2 r^^id 
$(P) oc 1/[(P/P*)^i + (P/P*)"'^] respectively. 

Our results are shown in Table |31 which gives the sam- 
ple in column (1), the evolutionary model in column (2), 
the two slopes (if applicable) of the LF in columns (3) and 
(4), the best- fit evolutionary parameter in column (5), 
and the break power (if applicable) in column (6). Errors 
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Fig. 8. — Maximum likelihood confidence contours (Icr, 2cr, and 
3(t) for the evolutionary parameter ki^ and the single power-law 
slope of the LF. The best fit values for the various classes are 
indicated by the various symbols, with the solid ones referring to 
the maximum likelihood analysis and the open ones pertaining to 
the Vo/Va approach. Note that the slope of the LF in the latter 
case depends on the adopted bin size in logarithmic power and is 
therefore only indicative. 

are la for one interesting parameter. The best-fit evolu- 
tionary parameters agree extremely well (mostly within 
la) with those derived through the K/K, approach. 

4. LUMINOSITY FUNCTIONS 

Fig. [5] shows that, assuming a single power- law LF 
and applying the maximum likelihood method, SFG and 
AGN, when considered as a single class, have widely dif- 
ferent forms of evolution, as seen above, but also different 
LF slopes, with $(P) oc p-^-^ and $(P) oc P'^-^ respec- 
tively. The 3(7 confidence contours have no overlap. The 
situation is even more extreme when one compares SFG 
with radio-loud AGN. On the contrary, in the case of the 
radio-quiet AGN the best fit parameters are the same 
as SFG (within less than Icr) and the confidence con- 
tours overlap to a large extent (although the radio-quiet 
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Fig. 9.— The local differential 1.4 GHz LFs for SPG and AGN 
in a P X 4'{P) form obtained with the maximum likelihood method 
(solid lines) and the LFs de-evolved to z = using the best 
fit evolutionary parameters from the maximum likelihood analy- 
sis (points). Filled circles indicate SFG, open squares radio-loud 
AGN, filled squares radio-loud AGN with P < lO^''-^ W Hz"!, 
and open circles radio-quiet AGN. The LFs for radio-loud AGN 
with P < lO'^*'^ W Hz~^ from both methods are shown without 
any de-evolution as inferred fro m the (Ve/Vk) analysis. Error bars 
correspond to la Poisson errors IIGehrelslll986l ) evaluated using the 
number of sources per bin with redshift determination only. See 
text for details. 

AGN contours are wider due to their smaller number). 
Finally, radio-quiet and radio-loud AGN have very differ- 
ent LFs and evolution parameters, with only a tiny over- 
lap of the 3cr confidence contours and the former having 
a much steeper LF than the latter ($(P) oc p-2.6 ^g^ 
$(P) oc p-i-45). 

The shape of the SFG LF is more complex than for the 
other classes, as the maximum likelihood double power- 
law fit, given in Table [Sj excludes the case 71 = 72 with 
a significance well above 3cr, which implies that a single 
power-law is not a good representation of the data. The 
best fit has $(P) ex P-^-^ and $(F) oc p-sis at the faint 
and bright end respectively, with a break at P '-^ 7 x 10^^ 
W Hz-i. 

Fig. ini shows the local LFs derived from the the max- 
imum likelihood best fits compared with those obtained 
from the l/Knax (hi^ our case 1/T4) technique without 
any assumption on the LF but de-evolved to z = us- 
ing the best fit evolutionary parameters. The LFs are 
shown in a P X $(P) form, which is almost equivalenlO 
to the 0(Mb) form normally used in the optical band 
and allows an easy separation of luminosity and density 
evolution as the former would simply translate the LF 



3 P X $(P) = 2.5//n(10) X $(A/) - 
the units of <J>(M) are mag^^ Volume^^. 



1.09 X *(M), where 
Note that the s e unit s 
are also sometimes used in the radio band: e.g., ICondonl 1)19891 ): 
ISadler et all pOOl ) : iMaiich fc Sadler! p007 '). 



to the right (higher powers) with no change in the ordi- 
nate (number), while the opposite would be true for the 
latter. 

The maximum likelihood fits, although relatively sim- 
ple, appear to be very good representations of the LFs 
obtained with the 1/14 technique. From Fig. [9] one 
also infers that AGN dominate over SFG for P > 3 x 
10^2 ^ Hz~^, in agre ement with previous studies (e.g., 
iMauch fc Sadleiil20Q7[ ). Moreover, radio-loud AGN have 
a much flatter LF than radio-quiet ones and are predom- 
inant at P > 6 X 10^^ W Hz^^. Finally, the radio-quiet 
AGN LF seems to be an extension of the SFG one at 
higher radio powers. 

It is important to be aware of the fact that the low- 
redshift behavior of the high-power radio-loud AGN is 
not well determined. Based on previous results (see Sec- 
tion [5321) '^6 would in fact expect a strong positive evo- 
lution at moderately low redshifts followed by a decline 
at higher redshifts. Our sample is too small to detect 
such a change, since it is dominated by high redshift ob- 
jects (~ 3/4 of the sample has z > 1). This means that 
the best fit PDF reflects largely the high redshift nega- 
tive evolution and, once the LF is de-evolved to z = 0, 
this translates into an artificially large density of sources 
at high powers. We therefore plot in Fig. [S]also the LF of 
P < lO'^'''^ W Hz^^ sources with their evolution fixed to 
zero, based on the results of the (14/14) analysis, which 
should give a more realistic estimate of the local LF. 

We now concentrate on the details of the individual 
classes. 

4.1. Star-forming galaxies 

Fig. [TOl shows different estimates of the local LF for our 
SFG. Filled triangles indicate the VLA-CDFS SFG LF 
in the 0.1 — 0. 35 redshift range, to compare it with the 
COSMOS LF (iSmolcic et al.l l2009aD in the same range 
(open triangles). The two LFs are within la apart from 
one bin. Our LF is "noisier" as we only have 24 objects 
in this redshift bin, as compared to 98 for the COSMOS 
sample. We also note that the selection criteria of the 
two samples are very different, with the COSMOS one 
being based on a re st-frame optical color classification 
(|Smolcicet al.ll2008[) . 

To have a more robust estimate of the local LF for 
SFG we have derived the LF de-evolved to z = (filled 
circles) using the best fit evolutionary parameter from 
the (14/14) analysis. This makes use of the whole sam- 
ple but is obviously dependent on the assumed evolu- 
tionary model. In Fig. [10] w e also plot th e best fits 
to the local (z < 0.3) L Fs from ISadler et al.l |2002) and 
IMauch fc Sadlej (|2007l ) (dotted and dash-dotted fines re- 
spectively) . Both our local LFs are consistent with these 
two LFs, particularly the de-evolved one {xt < 0.6). 
This validates our selection method. As was the case for 
the maximum likelihood approach, a single power-law fit 
is inconsistent with the de-evolved local LF {xf, ^ 3.4, 
significant at the ~ 97% level). 

The maximum likelihood fit provides a very good rep- 
resentation of the redshift evolution of the LF of SFG, 
as shown in Fig. [TT] (tabulated in Tab. [4]), which plots 
the SFG LF over the full redshift range sampled, that 
is — 2.3, in four redshift bins each containing a similar 
number of objects. The median enclosed volumes for the 
four bins are 3 x 10-'5,2 x 10-^9 x 10"^ and 5 x lO"'' 
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Fig. 10.— The local diflFerential 1.4 GHz LF for SFG in a PX(f>{P) 
form. Filled triangles indicate the VLA-CDFS LF in the 0.1 - 0.35 
redshift range, open triangles refer to the COSMOS LF in the same 
range, and filled circles denote the LF de-evolved to z = using 
the best fit evolutionary parameter from the (Vc/Vk) analysis. The 
best fits to the local SFG LF from Sadler et al. (2002) (converted 
to our value of Hq) and |Mauch & Sadler (2007) are also shown 
(dotted and dash-dotted lines respectively) . Error bars correspond 
to Ict Poisson errors (Gehrcls 1986) evaluated using the number 
of sources per bin with redshift determination only. See text for 
details. 



TABLE 4 
Luminosity Functions for VLA-CDFS star-forming 

GALAXIES 



Redshift range 


log Pi. 4 GHz 
WHz-1 


log P*(P) 
Gpc-3 




N 


0.038 <z< 0.214 


21.27 
21.67 
22.07 

22.47 


6.22l»;« 

t: cr;+0.25 
j.ao_g.2g 
t: CI +0.22 
O.Ot_0.24 






2 
9 
4 
5 


0.214 < z< 0.524 


22.19 
22.59 
22.99 


5 84+"'^^ 
t- OO+0.16 
^■'^'^-0.17 
4 40+0.29 






4 
9 
3 


0.524 < z< 1.021 


22.91 
23.31 
23.71 


c c:4+0.iV 
0.04-0.18 

4 Sfi+''-22 
4-00-0.24 

4.03l«il 


(5.44tj;;^34) 

(4.22«;fg) 
(3.90«:1) 


8(3) 
5(4) 
5(4) 


1.021 < z< 2.325 


23.79 
24.19 
24.59 


5.03±i!;i^ 

o 7n+0.20 
''■ '"-0.22 

^■'^-0.77 






10 
6 
1 


Numbers in parenthesis exclude the 
scale concentrations at 0.664 < z < 
l|Gilli ct al. 2003). Errors correspond 


seven sources in the two largc- 
0.685 and 0.725 < z < 0.742 
to la Poisson errors dGehreld 



[1986.) evaluated using the number of sources per bin with redshift 
determination only. The conversion to units of Mpc~ dex~ used, 
for example, bv ISmolcic et al.l (I2009ar . is done by subtracting 9 — 
log(ln(10)) from our values. 
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Fig. 11.— The differential 1.4 GHz LF for VLA-CDFS SFG in 
a P X (/'(P) form in four redshift bins. The solid lines represent 
the best fit double power-law LF from the maximum likelihood 
method evolved to the central redshift of the bin using the best-fit 
evolution (1 + z)^-*^, with dotted lines showing the same LF at 
the two extreme redshifts defining the bin. The short-dashed line 
represe nts the LF at z = 0. Error bars correspond to la Poisson 
errors (I Gehrelsl 119861) evaluated using the number of sources per 
bin with redshift determination only. The percentage of redshifts 
estimated from the optical magnitude is also given for each bin. 
Open symbols in the 0.524 — 1.021 bin do not include the seven 
sources in the two larg e-scale concentra tions at 0.664 < z < 0.685 
and 0.725 < z< 0.742 (IGilli et al.ll2003l) . See text for details. 



Gpc^ respectively. 



4.2. AGN 



Fig. [T2l shows different estimates of the local LF for our 
AGN. Filled triangles indicate the VLA-CDFS AGN LF 
in the 0.1 — 0. 35 redshift range, to compare it with the 
COSMOS LF (ISmolcic et al.l I2009bn in the same range 
(open triangles). The two LFs are within « la. As was 
the case for SFG, our LF is "noisier" as we only have 15 
objects in this redshift bin. The selection criteria for the 
two samples are also quite different and q uasars (< 20% 
of the A GN sample) were not included by ISmolcic et al.l 
((2009bl) . We also show the full LF for all AGN (filled 
circles) assuming no evolution, as inferred from our anal- 
ysis. 

In Fig. [12] we also plot the best fits to the local (z < 
0.3) L Fs from .Sadler et al.. (2002.) and .Mauch fc Sadleri 
()2007() (dotted and dash-dotted lines respectively). Both 
of our LFs are inconsistent with these previous estimates 
{xl ^ 3-0, significant at the > 98.9% level) and about a 
factor ^2-4 higher This is also true of the COSMOS 
LF of ISmolcic etall (|2009bl ) (xl ~ 4.6, significant at the 
> 99.9% level), despite their claims to the contrary. Two 
effects are at play here. The first, and more important 
one, is that the AGN LF includes a sizable contribution 
from radio-quiet AGN, which were not present in signif- 
icant numbers in previous LFs as these were based on 
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TABLE 5 
Luminosity Functions for VLA-CDFS radio-loud AGN 



22 24 26 

log P^ (1.4 GHz) [W Hz"'] 

Fig. 12.— The local diflFerential 1.4 GHz LF for AGN in a 
P X (t>{P) form. Filled triangles indicate the VLA-CDFS LF in 
the 0.1 — 0.35 redshift range, open triangles refer to the COSMOS 
LF in the same range, filled circles denote the full AGN LF (shifted 
by 0.1 in log P for clarity), while filled squares are the radio-loud 
AGN with P < lO'^*'^ W Hz~^ (with open squares showing the 
effect of excluding the ten sources in the two large-sc ale concen- 
tratio ns at 0.664 < z < 0.685 and 0.725 < z < 0.742: IGilU et ahl 
H2003I )). Both of these LFs are shown without any de-evolution as 
inferred from the (Vc/Va) analysis. The best fits to the local AGN 
LF from Sadler et al. (2002) (converted to our value of Hq) and 
[Mauch & Sadler (2007) are also shown (dotted and dash-dotted 
line s respectively). Error bars correspond to la Poisson errors 
IIGehrelail986i ) evaluated using the number of sources per bin with 
redshift determination only. See text for details. 

the NVSS survey (S'min ^ 2.8 mJy), while radio-quiet 
AGN make up a non-negligible fraction of radio sources 
only below « 1 mJy (Fig. H]), and both included only 
non-stellar optical sources. Note that none of these ar- 
guments apply to star-forming galaxies, which look non- 
stellar and, despite their relatively low fraction at high 
radio flux densities, were selected for because of the low 
redshift cuts (AGN tend to have higher redshifts: see Fig. 
[e]) and, in the case of the iMauch fc Sadleij (|2OO70 sam- 
ple, the K-ba nd selection. (This bias is vividly illustrated 
by Fig. 5 of iMauch fc Sadleil I2007L where star-forming 
galaxies dominate below ^ 10 mJy, while Fig. |4] shows 
that in purely radio selected samples without any redshift 
cut this happens below ^0.1 mJy.) If one considers only 
radio-loud AGN with P < lO^^-S w Hz^i (filled squares), 
which are non- evolving, our LF is marginally consistent 
with that from lSadler et all (|2002l) {xl ~ 2.0, significant 
at the ^ 95% level). The second effect is related to cos- 
mic variance: the exclusion of the eleven sources in the 
two large-scale concentrations at 0.664 < z < 0.685 and 
0.725 < z < 0.742 (jGiUi et al. 2003) (open squares), in 
fact, reduces even further the discrepancy with both pre- 
viously determined local LFs, which are now consistent 
with ours (xt < 1-8, significant at the < 91% level). 
The maximum likelihood best fits for the two classes of 



Redshift range 


log Pi. 4 GHz 


log P*(P) 


N 




WHz-1 


Gf 


C-3 






21.74 


5.79t«:f4 




3 




22.24 


4.591^7? 
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22.74 
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4 22+"-22 
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4.201^2! 
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24.24 
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1.0J_(j 20 


(4.351":^,) 


8(4) 
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s.mtlil 


(3.341^7?) 


4(1) 


0.651 <z< 0.964 
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4.iilSi^ 


i^-sitUl) 


6(3) 




25.42 


0-1+0.29 

o-»t_0.34 
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23.85 
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^•-^■^-0.22 
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24.35 


n on + 0.22 

J-»J_0.24 
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0.964 <z< 1.546 


24.85 


3A9+°ll 




4 




25.35 


3.361^3^ 




3 




25.85 


2.891^7? 




1 




23.82 


-|7+0.b2 
•^•-■-'-0.77 




1 




24.32 


., r,-|+0.36 




3 




24.82 


■^■^'-0.28 
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1.546 <z< 5.818 


25.32 


1-981^7? 




1 




25.82 


r, co+0.52 




4 




26.32 


1.961^7? 




1 




26.82 


2.44j:°;f, 




3 




27.32 


i.mtlfr 




1 


Numbers in parenthesis exclude the 


eleven sources in the two 


large- 


scale concentrations 


at 0.664 < z < 


0.685 and 


0.725 < z < 


0.742 


dGilli et al .1120031). Errors correspond to la Poisson errors dGehreld 


ll986D evaluated usine the number o 


■ sources per bin with redshift 


determination only. 


The conversion 


to units of 


Mpc-=^ dex-i 


used. 


for example, by |Smolcic et al.| (|2009b|), is done 


by subtractin 


g 9- 



log(ln(10)) from our values. 

radio-loud AGN, below and above P = lO^^-S W Hz"!, 
provide a good representation of the evolution and LF 
of radio-loud AGN, as shown in Fig. [13] (tabulated in 
Tab. [5]), which plots the radio- loud AGN LF over the 
full redshift range sampled, that is 0.1 — 5.8, in four 
redshift bins each containing a similar number of ob- 
jects. (For the radio-loud AGN with powers below lO^'*'^ 
W Hz~^ we fixed the evolution to zero, based on the 
results of the Vc/Va analysis, excluded the sources in 
the two large-scale redshift concentrations, and obtained 
^(P) (X p-i-^io.i-) ^]-^g median enclosed volumes for 
the four bins are 2 x lO"'', 7 x 10"^ 2 x 10"^ and lO'^ 
Gpc^ respectively. 

We note that the high-power end of the radio lumi- 
nosity function is poorly sampled because of the rela- 
tively small volume covered by our small-area field. This, 
combined with the fact that ^ 30% of the redshifts for 
our high-power radio-loud AGN have been estimated 
from the optical magnitude, makes our results for the 
P > 10^"''^ W Hz~^ sub-class more uncertain than for 
the other classes. Moreover, as discussed above, their 
low redshift evolution is not well determined (the two 
lowest redshift bins include two and four objects respec- 
tively), which also means that the best- fit LF derived 
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Fig. 13.— The differential 1.4 GHz LF for VLA-CDFS radio- 
loud AGN in a P X 4>{P) form in four redshift bins. The solid 
lines represent the best fit single power-law LF from the maximum 
Ukelihood method for P < lO^^'^ W Hz"! (left) and P > lO^^-S W 
Hz~^ (right), the latter evolved to the central redshift of the bin 
using the best fit for pure density evolution {l + z)~^'^ , with dotted 
lines showing the same LF at the two extreme redshifts defining 
the bin. The short-dashed line represents the best fit LF at z = 
for P > IQ-^^'^ W H z~^. Error bars correspond to Itr Poisson 
errors l|Gehrel£lll98(D evaluated using the number of sources per 
bin with redshift determination only. The percentage of redshifts 
estimated from the optical magnitude is also given for each bin. 
Open symbols in the 0.651 — 0.964 bin do not include the eleven 
sources in the two largo-scale concentrations at 0.664 < z < 0.685 
and 0.725 <z< 0.742 (Gilli et al. 20C)3). See text for details. 

by the maximum likelihood method is artificially high at 

5. DISGUSSION 

5.1. Effect of missing redshifts on our results 

Since 14/16 of the sources without redshift informa- 
tion are classified as radio-loud AGN (mostly based on 
their q and R values), the other classes are basically un- 
affected by our redshift incompleteness. We checked how 
the redshifts estimated from the optical magnitude influ- 
ence our results in two different ways. First, we assumed 
z = (z) = 1.13 for these 16 sources. This is very dif- 
ferent from our previous assumption, which resulted in 
a very broad redshift distribution extending between 0.3 
and 5.8. The (Vc/Va) (and best fit evolutionary parame- 
ters) changed negligibly for all but one sub-class, at most 
by an amount equal to OAa. Given the smaller redshifts 
(and therefore luminosity) involved, in fact, the sample 
of radio-loud AGN with P > lO^^-S W Hz"! shrank by 
~ 18%, while their (14/T4) decreased even further by 
1.5a reaching ~ 0.21. Second, we estimated the missing 
redshifts using the scatter of the correlation between log z 
and Vmag (dashed lines in Fig. [IJ. The anti-evolution 
of radio-loud AGN was confirmed despite the substan- 
tial redshift variations (Pev > 97.2%), with changes in 



{Vo/Va) only up to 0.6cr. In the case of high-power radio- 
loud AGN the changes were somewhat larger, with an 
increase in (14/14) of 1.6cr when the upward scatter was 
applied (which however implies estimated redshifts up to 
^ 10). The evolution was nevertheless still strongly neg- 
ative (Pov = 99%). As for low-power radio- loud AGN, 
the changes in (14/14) were < 0.3cr with still no signifi- 
cant evidence for a departure from the non-evolutionary 
case (Pcv < 86%). 

In summary, our results are quite insensitive to the spe- 
cific redshift values for the fraction of the sample without 
redshift information. This is not surprising, as redshift 
affects 14/14 values much less than flux density and our 
redshift incompleteness is very small (^ 



5.2. Effect of large-scale structures on our results 

I Gilli et aD ([2003) have studied the large-scale struc- 
ture in the CDFS in the X-ray and near-IR bands and 
detected two concentrations of sources in the 0.664 < z < 
0.685 and 0.725 < z < 0.742 ranges. Given the small area 
of our survey one could worry that such redshift spikes 
might influence some of our results. Indeed, the redshift 
distributions shown in Fig. |6]peak in the 0.5 — 0.75 bin 
for most classes but the peaks become much less pro- 
nounced when these objects are excluded. There are in 
fact 18 sources in these two redshift bins (7 SFG and 
11 AGN), which make up ~ 17% and ~ 11% (taking 
into account the effect of the sky coverage) of all SFG 
and radio-loud AGN respectively (but 0% of radio-quiet 
AGN and only ^ 3% of high-power radio- loud AGN). 

To assess the maximum impact of these two structures 
on our results we studied the evolution of our sources by 
excluding all sources in these two redshift bins. The re- 
sulting (14/14) values and best-flt evolutionary parame- 
ters were within la from those derived from the full sam- 
ples for all classes, which shows that the effect of these 
large-scale structures on our results is minimal. These 
over-densities are obviously more noticeable when one 
studies the evolution of the LF with redshift (see, e.g., 
the bottom- left panel of Fig. [TT] and the top-right panel 
of Fig. [T3| but even then the revised LFs are within < la 
from the old ones. The exclusion of these sources has also 
some effect when we limit the SFG and radio-loud AGN 
samples to 2; < 1.3 (Section [5] 



5.3. Are there two classes of low radio power AGN? 

We have identifled two classes of low-power AGN: 
radio-quiet ones, deflned as spelled out in Section [21 
~ 94% of which turn out to have P, < 10^^ W Hz^i, 
and radio-loud ones, characterized by Pr ^ 3 x 10^^ 
W Hz~^. Both classes have also relatively low R val- 
ues, as implicit in our sel ection of radio- q uiet s ources 
and as shown in Fig. 3 of iPadovani et all ([20091 ) . One 
obvious question is then what the differences between 
these two classes are. The answer is: many. First, they 
have very different (P > 99.99%) distributions in IRAC 
flux rati os (Section [2l), according to a two-dimensional 
KS test ( Fasano fc FranceschinilTl987f ). with radio- loud 
AGN mostly towards the old stellar population locus 
and the majority of radio-quiet AGN populating the 
region where most unobscured AGN should be (note 
that none of the selection requirements for the radio- 
loud/radio-quiet distinction was based on that). Sec- 
ond, despite the similar radio powers, they have very 
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different redshift distributions (P > 99.9%) and (z), 
1.73 and 0.84 for radio-quiet and fow-power radio-loud 
AGN respectively. Third, they evolve very differently 
and have different LPs, with radio-quiet AGN character- 
ized by a steep LP and strong positive evolution while 
low-power radio-loud AGN display a much flatter LP 
and no evolution (and in any case (T4/K) < 0.5; see 
Section [3]) . Pourth, they have very different X-ray-to- 
radio luminosity distributions and ratios. To properly 
take into account the upper limits on X-ray power we 
used ASURV (Lavallcy ct al. 1992), the Survival Anal- 
ysis package which emplo ys th e routines describ ed in 
iFeigelson fc NelsonI ()1985[ ) and llsobe et al.l (|1986D and 
evaluates differences in distributions and mean values 
by dealing appropriately with non-detections. The dis- 
tributions for the two classes are significantly different 
{P > 99.95%), with radio-quiet AGN having an order of 
magnitude larger X-ray power for the same radio power. 
Pifth, for the sources for which we have a spectrum the 
optical classification is very diff erent: ^ 2/3 of the ra dio- 
quiet AGN with a spectrum in iSzokolv et al.l ()2004D are 
either broad-lined or high-excitation line objects, while 
only ^ 8% of the low-power radio- loud AGN are classified 
as such (one object with Pr ~ 1.5 x lO^'' W Hz'^, which 
is close to the dividing line with radio powerful AGN). 
Pinally, even though the definition of the two classes is 
based on their q values, the point remains that low radio 
power AGN span ^ 3 orders of magnitude in far-IR to 
radio fiux density ratios (see Pig. [T4| , meaning that their 
radio emission goes from being related to star formation 
to having a likely jet origin. Pig. [T3] shows also that this 
is not the case for high-power (P^ ^ 3 x lO^'' W Hz~^) 
AGN, since they basically all have q < 1.7 or upper limits 
above this value. 

5.4. Are our radio-quiet AGN and star-forming galaxies 
really different classes? 

We have also identified two classes of radio sources with 
q values suggesting a star-formation origin for their ra- 
dio emission. Both classes have low R and low P^ (SPG 
by definition). They also evolve similarly in the radio 
band, have a similar slope of the LP, and the same range 
of q. Are we really dealing with two different classes? 
Again, the answer is: yes. Pirst, they have very dif- 
ferent (P > 99.99%) distributions in IRAC flux ratios 
(Section ^, according to a two-dimensional KS test, 
even before the few outliers were removed (P ^ 99%), 
with radio-quiet AGN populating the region where AGN 
should be and SPG distributed in the region where PAH- 
and starlight-dominated sources are expected to be. Sec- 
ond, despite their similar radio powers, they have very 
different redshift distributions (P > 99.99%) and (z), 
1.73 and 0.90 for radio-quiet AGN and SPG respectively. 
Third, for the sources for which we have a spectrum, the 
optical classification is very diff erent: ^ 2/3 of the ra dio- 
quiet AGN with a spectrum in iSzokoly et al.l (|2004[ ) are 
either broad-lined or high-excitation line objects, while 
none of the SPG are classified as such. Pinally, even 
though the definition of the two classes is based also on 
their X-ray powers, the point remains that radio sources 
with q > 1.7 span almost six orders of magnitude in 
Lx (~ three orders of magnitude on either side of the 
10**^ ergs s~^ divide). This means that their X-ray emis- 
sion goes from being related to star formation to having 




Fig. 14. — Radio power versus q, the logarithm of the ratio be- 
tween far-IR and radio powers for AGN in our sample. The ver- 
tical dashed line indicates q = 1.7, the dividing value between 
radio-quiet and radio-loud AGN, while the horizontal dashed line 
(Pr = 10^*'^ W Hz~^) divides low- and high-power radio-loud 
AGN. See text for details. 

a clear AGN origin. In summary, although radio-quiet 
AGN and SPG radio have similar q values and other 
radio-related properties, in one case we are clearly deal- 
ing with an AGN while in the other we are not. 

5.5. The evolution of micro- Jy radio sources 

We now analyze in more detail the evolution of the 
various classes and compare it with previous results. 



5.5.1. Star-forming galaxies 

Our results on the evolution of SPG in the radio band, 
which we fitted with P(z) oc (1 -f- z)''^ , Hl = 2.89to;i^ 
(or fci = 2.5l[!j from the (K/Va) method) in the 
range < z < 2.3, agree ve ry well w i th pr evious de- 
terminations. For example, iHopkind ()2004f ) obtained 
fci = 2.7 ± 0.6 up to z = 2 (and constant thereafter), 
with a small (but not significant) component of density 
evolution $(z) ex (1 -I- 2)Oi5±o.60^ rp^ check if we could 
constrain a possible redshift peak in the evolution we 
tried a simple model o f the type Pjz ] = ( 1 -I- z)^^^^ , 
first suggested by Wall. Pope fc ScottI (|2008f ). which al- 
lows for a maximum in the luminosity evolution followed 
by a decline. We found no evidence for /3 being signif- 
icantly different from 0, but we should stress that only 
two of our SPG have z > 2. 

As regards other bands, iMagnelli et al.l ([2003) have 
modeled the evolution of infrared luminous star-forming 
galaxies as a pure luminosity evolution P{z) ex (1-1- 



^) 



3.6±0.4 



up to z -- 1.3, while iMagnefliet al.l (|20lH ) have 



derived P{z) ex (l-f z)i"±0-9 for 1.3 < z < 2.3, which 
suggests a slowing down of the evolution at z > 1.3 (in 
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both cases the evidence for density evolution is not sig- 
nificant). If we split our sample we derive k^ — 3.5_q'7 
{{VJV^) = 0.63 ± 0.04, Pev = 99.98%) for z < 1.3 
(fci = 3.1!tio excluding the two large-scale structures) 

and kL = l.eta? ((^c/Ki) = 0.66±0.09, Pcv - 54%) for 
1.3 < z < 2.3, which is consistent with the IR results. In 
summary, although we find no significant evidence for a 
slowing down of the evolution at higher redshifts, such 
an occurrence cannot be ruled out by our data. 

Our SFG have (iin) ~ 5 x lO^^L©, which is typical 
for luminous infrared galaxies (LIRGs; 10^^ < Lm < 
IO^^Lq), and extend over 8 x 10^ < Ljr < 2 x IO^^Lq 
(where Lir refers to the 8 — 1,000 /im range and we 
ha ve used the mean va lue of the iiR/Pi.4Gffz ratio given 
in lSargent et al.|[2010l to estimate it), thereby reaching 
well into the ultra-luminous infrared galaxies (ULIRGs) 
regime (Lir > lO^^i©). 

5.5.2. Radio-loud AGN 

We have found a significant difference between the evo- 
lutionary properties of low-power and high-power radio- 
loud AGN, with a dividing line between the two at 
P ~ 3 X 10^"* W Hz~^. Namely, while low-power sources 
do not evolve, high-power ones anti-evolve significantly 
(P ~ 99.9%), which indicates that they were either less 
numerous or less luminous in the past. This is exactly the 
opposite of what was found in samples at higher flux den- 
sities (> 1 Jy), where high-power sources exhibit a strong 
positive evolution, with low-power ones still not evoly - 
ing fe.g.. lUrrv fc Padovanil[T99a IJackson fc Walllll999l) . 
This difference is likely due to two factors: 1. the range 
of redshifts sampled in the two cases are very different. 
For example, while a 10^^ W Hz~^ source in a sample 
defined by SiAGHz,min > 1 Jy can be detected only up 
to z '^ 0.065, if SiAGHz,m,in > 100 /iJy the same source 
can be seen up to z ~ 4.1 {a^ — 0.7). And indeed, for 
our high-power radio- loud AGN (z) ~ 2 (Tab. [5]), which 
is wh ere these sources are sup posed to start their decline 
fe.g.. lWaddington et al.|[200ll ): 2. a class of sources with 
P > 3 X 10^^ W Hz~^, moderate luminosity evolution, 
and a cutoff redshift Zmax = 5.5, will reac h a limiting 
flux density /min ~ 70 /zJy ()Padovanill201lD . When the 
flux density limit Smin of a survey is comparable to, or 
even lower than, this value, the observer will start "run- 
ning out" of sources and a deficit at higher redshifts will 
be observed. Stated differently, since at a first order 
Vc/Va ^ {S/Smin)~^^^, in this case S/Smin can often be 
> 1, if not ^ 1, which translates into small values of 
Vc/Va.. The key assumption here is that there needs to 
be a redshift cutoff, as if sources were present at all red- 
shifts then /min would be tending to zero and no such 
effect would be present. 

Note that the radio power, which separates the 
non-evolving from the anti-evolving radio sources co- 
incid es with the m inimum power of FR II sources 
(e-g.-lUrry fc Padov ani 1995t lGendre. Best fc Wall 201(1 
IPadovanill2011f ) and radio-quasars (IPadovanJ 20111 ). We 
then identify our low-power radio-loud AGN with low- 
luminosity radio galaxies of the FR I type and the high- 
power ones with the FR Il-like, powerful radio sources, 
which dominate the bright (> 1 mJy) radio sky, and of 
whic h we are w i tnessi ng the de mise. 

■Sadler et al.l (|2007D and .Donoso et al.l (|2009[) have 



studied the radio evolution of luminous red galaxies up 
to z — 0.7 and 0.8, finding evidence of weak but sig- 
nificant positive evolution. Their samples have very lit- 
tle overlap with ours, as they include radio sources with 
P > lO^'^-^- 1024-3 W Hz-i and only seven (~ 8%) of our 
AGN in this power range have z < 0.8 (see Fig. [S]). Fur- 
thermore, both studies only deal with red galaxies, while 
we include all radio sources. Our LF in the 0.4 — 0.8 
redshift range is in any case fully consistent with theirs, 
although our uncertainties are obviously quite large, once 
we exclude four sources belonging to the two large scale 
structures discussed in Section r5.2l and one optically com- 
pact source (which would have not been included by ei- 
ther of the two studies). The small, positive evolution 
they detect even at the lowest powers, that is between 
1024-2 and 1025-8 W Hz-i and 102^-3 and 102^ W Hz"! 
respectively, is consistent with these ranges encompass- 
ing the minimum value for FR II radio galaxies, which 
are known to evolve positively at the radio flux densi- 
ties s ampled by both studies (e.g.. lGendre. Best fc Walll 

ITavlor et al.l (|2OO90 have found that the number den- 
sity of massive (M^, > lO^^M©) red galaxies declines with 
redshift as $(z) oc (1 + 2)-i-60±o.i4(±o.2i) ^^ ^ < ^^g^ 
This is tantalizingly similar to the dependence we find for 
our high-power radio-loud AGN $(z) oc (1 -I- ^)-i-8±o.6 
over a larger redshift range (0.5 — 5.8). If we restrict our 
sample to z < 1.8 we still find evidence of negative evo- 
lution but not significantly so due to the smaller sample 
size. 

5.5.3. Radio-quiet AGN 

We have estimated for the first time the evolution of 
radio-quiet AGN in the radio band, modeling it as a pure 
luminosity evolution and obtaining P{z) ex (1 -I- z)'^^, 
kL = 2.5to:5, in the range 0.2 < z < 3.9. In the X-ray 
band the situation appears to be more complex, with 
strong evolution up to z ~ 1 — 2 ( depending on lu minos- 
ity) and then a slow down (e.g., iHasinger et al.| [2005'). 
As done for SFG, we then tried a model of the type 
P{z) = (1 -I- z)*'+'^^ but our previous best fit {k = 2.5, 
13 = 0) was still fully consistent with the data, although 
peaks at z « 3, for example, are also within the la 
contours. Larger samples of radio-selected, radio-quiet 
samples will put strong constraints on the evolution of 
these sources, by-passing also the problems of obscura- 
tion, which plague the optical and soft X-ray bands, al- 
though source identification will require ancillary multi- 
wavelength data, as shown in Section [2j Moreover, the 
surface density of our radio-quiet AGN, ~ 520 deg~2^ 
is already a factor ^ 6 larger than that of one of the 
deepest opt ically-selected quasar samples (~ 80 deg~2 to 
g ~ 22; Ifilchards et al. 2005) and only - 1/4 of that of 
unabso rbed AGN down to fa int Chandra fiuxes (~ 2, 000 
deg~2; IHasinger et al.|[2005l) . (Note however that we are 
sensitive to both broad- and narrow-lined AGN, while 
both the optical and X-ray samples under consideration 
include only the former type.) 

Our radio number counts are consistent with those in 
the hard X-ray b and, as shown i n Fig . 31 which shows 
the predictions of iWilman et al.l (|2008[ ). based on a con- 
version of the AGN X-ray LF to a radio LF, and of Pa- 
per IV, obtained from the X-ray number counts by using 
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a typical radio-to-X-ray flux density ratio. This shows 
that the sources we are selecting in the radio band are 
the same as the X-ray emitting radio-quiet AGN. 

5.6. Comparison with the COSMOS and Deep SWIRE 
Field surveys 

iSmolcic et all (|2009al ) and lSmolcic et all (|2009bf ) have 
studied the evolution of SFG and AGN in the radio band 
using 1.4 GHz VLA observations of the COSMOS 2 deg^ 
field, which ha ve a limiting flux densi ty ~ 45 /zJy in the 
central 1 deg^. iStrazzullo et al.l (|2010[ ) have studied a 1.4 
GHz selected sample of radio sources in the Deep SWIRE 
Field (DSF), reaching a limiting flux density ~ 13.5 fiJy 
at the center of a 0.36 deg^ area. We can then compare 
our results more directly with theirs, keeping in mind 
that both samples reach only z = 1.3 and that quasars, 
which make up < 20% of the total, were not included in 
t he AGN COSM OS sample. 

ISmolcic et all (|2009a) derived fc^ = 2.1 ± 0.2 or fci, = 
2.5 ± 0.1 for SFG depending on the choice of the local 
LF, while we get /cl = S.bto^ for z < 1.3 (fc^ = 3.lt% 
excluding the two large-scale structures). Their values 
are smaller than ours but not signiflcantly so, given 
our relatively large uncertainties and the presence of 
large-scale structures in our field. Their evolution, un- 
like ourSj_Js_Jiowevei_^gnificantly weaker than found 
by iMagnelhetan (|2009f ) in t he IR b and for z < 1.3 
k r, = 3.6 ± 0.4: see Section I5.5.1|) . IStrazzullo et al.l 



20101) have defined two classes of "star-forming" (blue) 
and "intermediate" (green) galaxies, for which they de- 
rive fci = 2.9 ±0.3 and fci = 3.6 ±0.2 respectively (their 
non-param etric results). This evo lution is stronger than 
obtained bv lSmo lcic et al.' ('2009a') but in good agreement 
w ith ours and a lso MaancUi ct al. (2009). 

ISmolcic et alj (|2009b) found fci = 0.8 ± 0.1 or fci, = 
1.1 ±0.1 for their AGN. For z < 1.3 we obtain (K/K) = 
0.54 ±0.03, indicative of slightly positive but not signifi- 
cant (Pev = 77%) evolution, with ki, — 1.5 and kjj = 1.0 
(and obviously large error bars). If we exclude the two 
large-scale structures (K/K) = 0.57 ± 0.03, indicative 
of positive (Pcv = 96.9%) evolution, with k^ = 2.3 ± 1.0 
and kjj = 2.0 ± 1.0. These are larger than the COS- 
MOS values but still consistent with them given our rel- 
atively large error bars. Part of the difference might also 
be explained by the exclusion of quasars from the COS- 
MOS sample, as these are expected to be mostly of the 
radio-quite type (Paper IV) and therefore strongly evolv- 
ing. Indeed, when splitting the sample into radio-quiet 
and radio-loud AGN one obtains (Ve/K) = 0.69 ± 0.06, 
kL = 3.9tg-^ (Pev ~ 96.8%) and (K/K) = 0.49 ± 0.04 
(Pov ^ 85%) respectively ((K/K) = 0.52 ± 0.04 ex- 
cluding the two large-scale structures). The apparently 
positive AGN evolution appears then to be driven by 
the radio-quiet sources, while the radio-lo ud ones do not 
evolve in the redshift range sampled by ISmolcic et al.l 
(|2009bf) . This has important implications for their re- 
sults, as they do not distinguish between the two classes 
of AGN and use the positive evolution they found to 
estimate, for example, the evolution of the comoving 
radio luminosity density and mechanical energy output 
of low-radio-power AGN. IStrazzullo et al.l ()2010[ ) derived 
/cl = 2.5 ±0.3 for their "quiescent" (red) galaxies, which 
they identify as AGN, which is in good agreement with 



our result but significantly larger tha n the COSMOS 
value. Note that -- 40% of the AGN in iStrazzuUo et al.l 
(|2010f) are within la of the infrared-radio correlation typ- 
ical of SFG, which shows that, like ours, about half of 
their AGN cannot be of the radio-loud type (see also 
Prandoni 2010 ). 

5.7. The origin of radio emission in radio-quiet ACN 

The mechanism responsible for the radio emission 
in radio-quiet AGN has been a matter of debate 
since the discovery of quasars. Alternatives have 
included a scaled dow n version of the radio-loud 
AGN mecha nism (e.g.. iMiller et al.l Il993( ). star for- 
mation (e.g., iSopp fc AlexandeHll99lD . a magnetically 
heat ed corona ("Laor fc Behaii 120081). and disk winds 
(JBlundell fc Ku ncic 20Q3 (biat see iSteenbrugge et all 

Our results suggest very close ties between star forma- 
tion and radio emission in radio-quiet AGN, since their 
evolution is indistinguishable from that of SFG (Sec- 
tion [3]) and their LF appears to be an extension of the 
SFG LF (Section H]). Furthermore, radio emission in the 
two classes of AGN is bound to have a different origin. 
If radio-quiet AGN were simply "mini radio-loud" AGN, 
in fact, they would have to share the evolutionary prop- 
erties of the latter and their LF should also be on the 
extrapolation of the radio-loud one at low powers. None 
of these two facts is borne out by our data (see, e.g.. Fig. 

El). 

This concurs with the re sults of various papers 
over the past 20 y ears (e.g., iSopp fc Alexanden 119911 : 
I Sargent et al.l I2010D , which have shown that radio- 
quiet AGN and star-forming galaxies have very similar 
FIR- to-radio flux density ratios. Note, however, that 
Sand ers et al.l (l989) ha ve discarded this as e ntirely co- 
incidental, and similarly iKukula et al.l (|1998( ) have sug- 
gested that dust heating by the quasar and AGN-related 
radio emission could also conspire to make this hap- 
pen. Moreover, the detection of compact, high bright- 
ness temperature cores in several radio- quiet AGN (e.g., 
lUlvestad. Antonucci. fc Barvainisll2005[ ). which resemble 
those observed in radio-loud AGN, would also argue 
against our results, although in some Seyfert galaxies 
these cores are surrounded by d iffuse radio emission co n- 
ncctcd to star- forming regions (jOrienti fc Prietoll2010[ ). 

This suggests that AGN and star-formation related 
processes coexist in radio-quiet AGN. Indeed, the frac- 
tion of flux density contained in the compact, central 
component is, on average, ~ 70% for low-redshift radio- 
quiet AGN (e.g., iKellermann et al.lll989l: IKukula et all 
Il998f ). which leaves some room for extended emission. If 
the AGN-related component is non-evolving, as appears 
to be the case for low-power radio-loud AGN, while the 
star-formation-related one follows the evolution of SFG, 
one could understand the difference between our radio- 
quiet AGN, which have (z) ~ 1.7, and those imaged with 
the VLBI, which are mostly local. In fact, since this red- 
shift difference implies an increase by a factor « 10 — 20 
in the radio power related to star formation, a minor 
(^ 1/3) extended component at z ~ would then be- 
come dominant (> 3) at higher redshifts. Furthermore, it 
should be pointed out that no complete sample of radio- 
selected, radio-quiet AGN has ever been observed at the 
VLBI resolution. The choice of objects that could be 
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detected might have then led to a selection effect favor- 
ing relatively strong targets, more similar to radio-loud 
sourceS; 

Ric hards et al.l (|2007| ) have studied a sample of 92 ra- 
dio sources in the Hubble Deep Field North brighter than 
40 ^i3y and well resolved by MERLIN and the VLA. By 
classifying more than 70% of them as starbursts or AGN 
using radio morphologies, spectral indices, optical ap- 
pearance, and rest-frame MIR emission, they found that 
the X-ray luminosity indicates the presence of an AGN in 
at least half of the 45 radio starbursts with X-ray coun- 
terparts. Moreover, almost all extended radio starbursts 
at z > 1.3 host X-ray selected obscured AGN and their 
radio and X-ray powers are uncorrelated, which points to 
different emission mechanisms being at play in the two 
bands. These results, which associate high-redshift ra- 
dio starbursts with AGN, are fully consistent with our 
suggestion of a very close relationship between star for- 
mation and radio emission in our relatively high-redshift 
radio-quiet AGN. 

If radio emission in radio-quiet AGN is mostly related 
to star formation processes, we ca n use the mean value 
of the Liji/PiAGHz ratio given in I Sargent et al.l (|2010f ) 
to estimate their IR powers. Our radio-quiet AGN have 
(Lir) '--^ 4 X 1O^^L0, that is well in the ULIRGs regime, 
and reach 2 x lO^^'i©. The star- formation rates (SFRs) 
implied by their radio powers assuming that all radio 
emission is star-form ation r elated and using the relation- 
ship derived bv iBelll (|2003[ ). are, on average, ~ 500 Af© 
yr~^, extending over the '^ 10 — 20, 000 M© yr~^ range. 
The mean value is typical of ULIRGs, while the up- 
per end is in the hyperluminous infrared galaxies regime 
(LiR > 10 "L©), which can reach SFRs > 10,000 M© 
yr~^ (e.g.. lRowan-Robinson fc Wandl2010D . 

The association of ULIRGs with radio-quiet AGN is 
certainly not new. I Sanders et al.l fl988) proposed an 
evolutionary connection between ULIRGs and quasars, 
based on the fact that all of the twelve ULIRGs in 
their sample displayed AGN spectra in the optical band 
(~ 2/3 of their quasars are radio-quiet). Moreover, the 
AGN detection rate amongst local ULIRGs is ^ 70% 
([Nardini et al.l I2010D and radio-detected ULIRGs are 
known to be rare locally but rap idly evolving with red- 
shift re.g.. lMauch fc Sadlei<l2007h . 

5.8. Astrophysics of micro- Jy sources 

'Fanaro ff fc Rilevl (|1974[ ) recognized that radio galaxies 
separate into two distinct luminosity classes, each with 
its own characteristic radio morphology. High-luminosity 
FR lis have radio lobes with prominent hot spots and 
bright outer edges, while in low-luminosity FR Is radio 
emission is more diffuse. The luminosity distinction is 
fairly sharp at 178 MHz, with FR Is and FR lis ly- 
ing below and above, respectively, the fiducial luminos- 
ity Pi78MHz ~ 1026/(ffo/70)2 W Hz"!. This translates 
at higher frequencies to Fi.4ghz « 3 x 10'^^'/{Ho/70)^ W 
Hz""'^ (assuming a^ = 0.7), w ith some dependency also on 
optical luminosity (jOwen fc W hite 1991) and therefore a 
rather large overlap. An independent separation on the 
basis of nuclear activity into high-excitation (HERGs) 
and low-excitation (L ERGs) radio-galaxi es has been pro- 
posed more recently (iLaing et al.lll99l) . It turns out 
that almost all FR Is are LERGs and most FR lis are 



HERGs, although there is a population of FR II LERGs 
as well. Observational evidence indicates that the two 
types of radio galaxies have intrinsically different cen- 
tral engines. Namely in LERGs the accretion disk, if at 
all present, is thought to be m uch less efficient than in 
HERGs (e.g iChiaberge et al.| [T999: Donato et al. 200i 
lEvans et al.ll2006D . This points to a large difference in 
a ccretion rates b etwee n the two classes. 

iCroton et al.l ()2006[) have associated high-accretion 
sources, and therefore also HERGs, with their so-called 
"quasar-mode" , which they interpret as merger-driven, 
efficient accretion of cold disk gas, present also in radio- 
quiet AGN. Low- accretion objects, i.e., LERGs, on the 
other hand, have been connected with the less efficient 
accretion of warm gas, the so-called "radio-mode" . The 
mean black hole accretion rate for the "radio-mode" is 
predicted to be approximately constant up to 2: ss 2, 
based on a suite of semi-analytic models implemented 
on the output of the Millennium Run (Fig. 3 of 
ICroton et al.ll2006l ). "Quasar-mode" accretion, on the 
other hand, is envisioned to be most efficient at z ^ 2 — 4, 
dropping by a factor of 5 by z ~ 0. This is similar in form 
but somewhat weaker than the observed cosmological lu- 
minosity evolution of br ight quasars in th e radio, optical, 
and X-ray bands (e.g., iWall et all 120051 and references 
therein) . 

As we identify our low-power radio-loud AGN with FR 
Is/LERGs these would then be inefficient accretors in a 
"radio-mode" , which would explain their lack of evolu- 
tion, assuming that it is driven only by the accretion 
rate. We instead identify as sources in a "quasar-mode" 
the high-power radio-loud AGN (which should be mostly 
HERGs) and the radio-quiet ones. Their different evolu- 
tions in the radio band are explained by distinct emission 
mechanisms and by the fact that the high radio powers 
of the former mean that we are seeing the effects of a 
high redshift cutoff. 

6. SUMMARY AND CONCLUSIONS 

We have used a deep, complete radio sample of 193 ob- 
jects down to a 1.4 GHz flux density of 43 /iJy selected in 
the Chandra Deep Field South area to sharpen our un- 
derstanding of the nature of sub-millijansky sources and 
to study for the first time their evolution and luminosity 
functions up to z ~ 5. Our unique set of ancillary data, 
which includes far-IR, near-IR, and optical observations, 
redshift information, and X-ray detections or upper lim- 
its for a large fraction of our sources, has allowed us to 
develop an unprecedented classification scheme to cate- 
gorize in a robust way faint radio sources in star-forming 
galaxies, radio-quiet, and radio-loud AGN. Our main re- 
sults can be summarized as follows; 

1. Star- forming galaxies and AGN make up an equal 
part of the sub-millijansky sky down to 43 /xJy, 
with the former becoming the dominant population 
only below w 0.1 mJy. Radio-quiet AGN are con- 
firmed to be an important class of sub-millijansky 
sources, accounting for ~ 30% of the sample and 
~ 60% of all AGN, and outnumbering radio-loud 
AGN at < 0.1 mJy. 

2. The radio power of star- forming galaxies evolves 
as (1 -|- z)'^-^~'^-^ up to z < 2.3, their maximum 
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redshift in our sample, in agreement with previ- 
ous determinations in the radio and IR bands. Al- 
though evidence of a slowing down of the evolution 
at z > 1.3 is not significant it cannot also be ruled 
out. The radio luminosity function of SFG can be 
parametrized as broken power-law $(P) oc P~^'^ 
and ^{P) oc p^3i5 at the faint and bright end re- 
spectively, with a break at P ~ 7 x 10^^ W Hz~^, 
which is also consistent with previous derivations. 

3. AGN as a whole do not appear to evolve. However, 
once they are split into radio-quiet (energy bud- 
get dominated by thermal emission) and radio-loud 
(dominated by non-thermal, jet emission) sources 
the situation is very different, with the radio-quiet 
population evolving very significantly and simi- 
larly to star-forming galaxies and the radio-loud 
population displaying negative density evolution 
^(z) oc (1 + z)~^-®='=°-^. The luminosity function of 
radio-loud AGN, $(P) oc P"^-^, is also much flat- 
ter than that of radio-quiet AGN, $(P) oc P'^*^, 
derived here for the first time, which seems to be 
an extension of that of star-forming galaxies at the 
high power end. 

4. There is a significant difference between the evolu- 
tionary properties of low-luminosity radio galaxies 
and radio powerful (F > 3 x 10^^ W Hz"!) AGN, 
as while the former do not evolve, the latter evolve 
negatively. This is exactly the opposite of what 
found in samples at higher flux densities (> 1 Jy), 
where high-power sources exhibit a strong positive 
evolution. We interpret this difference as due to 
the fact that we are sampling the high-power radio- 
loud population up to much larger redshifts (z > 5) 
and as a result of a redshift cutoff. 

5. Our results suggest a very close relationship be- 
tween star formation and radio emission in radio- 
quiet AGN, since their evolution and luminosity 
function are respectively indistinguishable from, 
and an extension of, that of star-forming galax- 
ies. This is supported by the fact that radio-quiet 
AGN and star-forming galaxies appear both to fol- 
low the "IR-radio relation" but is in contradic- 
tion with the detection of compact, high brightness 
temperature cores in several (mostly local) radio- 
quiet AGN, similar to those observed in radio-loud 
ones. The co-existence of two components, one 
non-evolving and AGN-related, and one evolving 
and star-formation-related, and selection effects in 
the choice of radio-quiet VLBI targets, can recon- 
cile these apparently discrepant observational data. 

6. The surface density of radio-selected, radio-quiet 



AGN, ~ 520 deg~^, is already about 6 times larger 
than that of one of the deepest optically-selected 
quasar sample and only ~ 1/4 that of unab- 
sorbed X-ray selected AGN. This means that sub- 
millijansky radio surveys, given the appropriate an- 
cillary multi-wavelength data, have the potential 
of detecting large numbers of radio-quiet AGN by- 
passing the problems of obscuration, which plague 
the optical and soft X-ray bands. The radio num- 
ber counts of radio-quiet AGN are consistent with 
those in the hard X-ray band, which shows that the 
sources we are selecting are the same as the X-ray 
emitting radio-quiet AGN. 

7. Sub-millijansky radio surveys wanting to study 
the evolution of faint, radio-loud AGN, need to 
consider that a large fraction {^ 60% down to 
~ 50 /iJy) of the radio-selected AGN are actually 
of the radio-quiet type and therefore need to be 
treated separately. This has not been done so far 
by other studies and can have a large impact on the 
study of faint "radio-mode" inefficient accretors. 

We plan to expand on this work by using our deeper 
radio observations (M iller et al.l 120081 and in prepara- 
tion) and the re cently released 4 Msec Chandra data 
()Xue et al.ll201lD . This will provide us with a catalogue 
of ~ 900 radio sources, with which we will be able to ad- 
dress the issues discussed in this paper with larger statis- 
tics. 
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